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1 Introduction 



In this paper we derive and give examples of a generalized Fourier series for complex binomials of the form [z — cos V'l ^^i 
where ^SR, /iSC, 2:SC \ (—00, 1] and \z\ > 1. We show that the Fourier series is given compactly with coefficients 
given in terms of odd-half-integer degree, general complex-order associated Legendre functions of the second kind. 
(Hereafter we refer to associated Legendre functions simply as Legendre functions.) Algebraic functions such as 
this naturally arise in classical physics through the theory of fundamental solutions of Laplace's equation, and they 
represent powers of distances between two points in a Euclidean geometry. 

Fourier expansions for algebraic distance functions have a rich history, and this expansion makes its appearance 
in the theory of arbitrarily-shaped charge distributions in electrostatics ([33], [H], [S], [35], [31], [10] )j magnetostatics 
([35], [B]) quantum direct and exchange Coulomb interactions ([H], [10], [IS], [32], [1]), Newtonian gravity ([H], [37] . 
[26] . [25] . [I9] . [31], [36], [9], [28], [7], [38]), the Laplace coefRcients for planetary disturbing function ([l4]. [I5]). and 
potential fluid flow around actuator discs ([8], [24]), just to name a few direct physical applications. A precise Fourier 
^im4> analysis for these applications is extremely useful to fully describe the general non-axisymmetric nature of these 
problems. 

The fact that this classical theory becomes recently relevant is due to several reasons. Firstly the fundamental 
solutions for Laplace's equation are ubiquitous in Pure and Applied Mathematics, Physics, and Engineering. Secondly 
the fact that Fourier expansions encapsulate rotationally-invariant symmetries of geometrical shapes makes it an ideal 
model case for critical study above and beyond the purely spherically symmetric shape, and therefore it provides 
a powerful tool when implemented numerically in a variety of important problems. Thirdly, Fourier expansions for 
fundamental solutions gives rise to a hardly-studied aspect of Special Function Theory, which allows one to further 
explore the properties of higher transcendental functions. 



2 Generalized Heine identity 

Gauss (1812) ([2T], p. 128 in Werke III) was able to write down closed-form expressions for the Fourier series for the 
related function [rf -l-r| — 2rir2 cos ■(/']"''. Gauss recognised that the coefficients of the expansion are given in terms of 
the 2F1 hypergeometric function, and he was able to write down a closed-form solution (where we have used modern 
notations) given by 
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[rl + rl — 2rir2 cosi/']^ 



^.„cos(nV) — ■ 271+;;- 2Fi[n + fi,^;n + l;-j], (1) 



assuming ri,r2 G M, r2 < ri, e„ is the Neumann factor [29| e„ = 2 — (5„ where (5„ is the Kronecker delta, commonly 
occurring in Fourier expansions, (/i)„ is a Pochhammer symbol for a rising factorial (see eqs. (O and (O), and the 
Gauss hypergeometric function is defined through 



(«)fe(&) 



k=0 ^ 



Carl Neumann (1864) ([30], P- 33-34 therein) was one of the first to study separable solutions to Laplace's equation in 
toroidal coordinates. He examined and wrote down the Fourier expansion in term of Gauss hypergeometric functions 
for the instance /i = i. 

But it wasn't until Heine (1881) in his Handbuch der Kugelfunctionen ([23], p. 286 therein) that it was recog- 
nised that this particular Gauss hypergeometric function represented a certain special class of higher transcendental 
functions, namely Legendre functions of the second kind with odd-half-integer degree. These Legendre functions, 
and in particular, those with integer order, are toroidal harmonics, the functions which separate Laplace's equation in 
toroidal coordinates. The fact that the algebraic function of present study relates to behaviour on toroids of revolution 
is an important consequence of a non-axisymmetric Fourier description of potential theory in rotationally-invariant 
coordinate systems. We now proceed to build upon Heine's formula in order to derive a generalization of his identity. 

2.1 Derivation of the identity 

We are interested in computing the following Fourier expansion 
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1 °° 

— = ^cos(n7A)A^,„(z), (2) 



[z — COS V'] 

where G C, z G C \ (—00, 1], and \z\ > 1. 

The expression for these Fourier coefRcients is given in the standard manner by 

A f \ ^" r C0S("V') ^ , 



COS ipY 

We can collapse to purely non-negative modes if the summand is an even function of the Fourier quantum number. 
In that case 

00 00 
^ ^^,„(z)e™'^ = ^ enAf,,^{z) cos(nV'). 



71— — 00 n— 

In order to do this, we must show that 

/ cost/iX 

arg (z — cos ip) — arg(z) + arg II 1 . (3) 

Eq. ([3]) is verified as follows. We define the function / : [— 1, 1] ^ M such that 

f{x) = arg {z + x)~ arg(z) - arg (l + ^) • 

/ is clearly continuous and since [—1,1] is connected, /([— 1,1]) must be connected and /([— 1,1]) C 27rZ. Hence 
/([—1, 1]) is a one-point set and since /(O) = 0, / is a constant equal to zero. Therefore we have shown that eq. ([3]) is 
true and hence we can rewrite the left-hand side of eq. ^ without loss of generality as 



[z-costpf 



cos-^ 



Since \z\ > 1 and cosi/; S [—1,1] this implies that \{costp)/z\ < 1, and we are in a position to employ Newton's 
binomial series 

where G C, \w\ < 1, and 

V\ 1 r(/i + 1) 



^kj k\T{fi-k + l) 

is the generalized binomial coefficient. 

In what follows we utilize Pochhammer symbols representing either rising or falling factorials. Unfortunately 
there is no standard convention used for Pochhammer symbols. So here we utilize the following convention which is 
consistent with usage in Special Function theory (i.e. with hypergeometric functions). The Pochhammer symbol for 
rising factorial is given by 

f 1, if n = 0; 

iz)n^l (5) 

\ (z)-(z + l)-(z-f2)---(z + n-l), if n>l, 

where n € Nq = N U {0} and z G C. Note that the Pochhammer symbol for the rising factorial is expressible in terms 
of Gamma functions as 

Similarly the Pochhammer symbol for the falling factorial, with the same quantities, is given by 

[1, if n = 0; 

\z\n — \ 

1 (z) • (z- 1) • (z-2)---(z-n-f 1), if n>l. 
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The Pochhammer symbol for the faUing factorial is also expressible in terms of Gamma functions as 

^ ^" " T{z-n + iy 

One can clearly relate the two symbols as such 

(-l)"hz]„ = (z)„. 

The binomial coefficients can be defined through the Pochhammer symbol for the falling factorial 

z\ _ [z\n 
nj n\ 

Combining this relation for the binomial coefficients, eq. ([7]), and eq. (|4]), we obtain 



(7) 
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[z — cost/i]^ fc! 



E 



-fi — k 



COS Ip. 



(8) 



We can expand the powers of cosine using the following trigonometric identity 

k 

2k 



1 



n=0 



cos[(2n — k)tp], 



which is the standard expansion for powers using Chebyshev polynomials (see for instance p. 52 in Fox and Parker 
(1968) [18|V Inserting this expression in eq. ([5]), we obtain the following double summation expression 



oo k 

EE 



^ ^' k=On=0 



cos[(2n — k)tp]. 



(9) 



Now we perform a double index replacement in eq. There are two separate cases k <2n and k > 2n. There is 
an overlap if both have an equality and in that case we must multiply by 1/2 after we sum over both cases. If A: < 2n 
make the substitution k' = k — n and n' = 2n — k. It follows that k = 2k' + n' and n = + fc', therefore 



2fc' + n' 
n' + k' 



2k' + n' 
fc' 



If fc > 2ri make the substitution k' = n and n' = k — 2n. Then fc = 2fc' 



and n ~ k' therefore 



2fc' + n' 
k' 



2k' + n' 



where the equalities of the binomial coefhcients are confirmed using the identity 



I 

I — m J ^ 



(10) 



where I d Ij. To take into account the double counting which occurs when k = 2n (which occurs when n' = 0) we 
introduce a factor of e„'/2 into the expression (and relabel k' ^ k and n' n). We are left with 



cos ■0]^ 



2zf' 



^ OO 

— ^ e„ cos(nV) 

(M)2fe+n 



n=0 

oo 

E 

fc=0 



(2fc + n)! (2z)2fe+'^ 



2fc + n 
fc 



2fc + n 
k + n 



which is straightforwardly simplified using the definition of the binomial coefhcients and eq. (|T0 
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cos iJjY 



-j^ oo oc 

— > en cos(n^h) ——— y 



2k+n 



1/1 



ri=0 



(22)"^fc!(fc + n)!4'= 
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The second sum is given in terms of a Gauss hypergeometric series 



- 5: e„ cosM)-^i^ ( ; n + 1; - ) . (11) 



r!=0 ^ 



This is just a re-statement of Gauss's original formula, eq. ([1}, with 

only now the quantities z, /i G C. This Gauss hypergeometric function is expressible in terms of the Legendre function 
of the second kind (Chapter 8 by I. A. Stegun [1], p. 332) 

2F1 I — : — , z ; V + 1: 



2 2-T{u + I)z'-+A'(z2 - l)-M/2+l/4g-».(p-l/2) 

where i^, /i € C, a branch cut is chosen along the real axis along (—00, 1], and otherwise this expression is valid except 
where the constituent functions are not defined. 

If we substitute = n g Z in the hypergeometric function and take advantage of the following property of Legendre 
functions of the second kind with |z| > 1 |13j 

= *3n-l/2(^)' 

in eq. (jlip . for all 71 G Z and G C, we obtain a complex generalization of Heine's reciprocal square root identity 
given as follows 

1 /2"e-'"'(^-l/2),' 2 _ 1^-p/2+l/4 °o 

= v!' hr ^ ^-or;;?(=). m 

where Heine's original identity ([13], p. 286) is the case ioi — \ given by 

1 \/2 



\/ Z — cos TT 



^ Q„_i(z)e™'^, (14) 



m— — 00 



see [H] for exact forms and recurrence relations for these Legendre functions. Eq. (|14p was recently generalized for /i 
given by odd-half-integers in Selvaggi et al. (2008) ([39] 7 P- 033913-6). 

The generalization given by eq. (jl3p is also expressible in terms of Legendre functions of the first kind through 
Whipple's transformation of Legendre functions ([13], [T]) as ,, 

1 ( 2 _ l\-p/2 °° 



T J \/ y n — — 00 ^ ^ 



In the following section we describe some specific examples of the generalized identity, and present some interesting 
implications. 

2.2 Examples and implications 

We now have the value of the following definite integral 

r cos{nt)dt „ ^ e-"(^-V^)(z^-l)-^/^+V^ ^^_i/,, 
[z-coBtY ~ ^ fV) y«-i/2(^)- 
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The principal example for the generalized Heine identity which was first proven in Selvaggi et al. (2008) [SS] (see 
also [17] I P- 182 and [W) for fi~ q + ^ where q e Z is given by 

[z-cos^]9+i/2 " 7r(2g - l)!!(z2 _ i)«/2 ^" cos(nV)Q«_i/2(z), 

where 2;!! is the double factorial (Chapter 6 by P. J. Davis 1 ). Notice that the double factorial is well-defined for 
negative odd integers 

^ ^ (2n)! ■ 

For instance, for q = —1 we have 



, - 1 En COs(n?A) 

- COS^ = ^ 2^ Qn-lM^), 
^ n— 4 



and for q = +1 

_23/2 

[z - cosV']^'^^ = 



= ^e„cos(nV')Q,\_i/2(2). (15) 



TTy/ — ^ „ 
^ n— 

Note that the minus sign in eq. (jlSp expansion might seem initially troublesome, except that it is important to notice 
that the unit-order Legendre functions of the second kind are all negative in sign as can be easily seen from the Gauss 
hypergeometric function representation, eq. (|12p . We have as such 



'/'^ ^ ^2(z-i) \y z + 1 



and 



'^'^ ^ ^2(z-i) [yz+ij V 2 \vz+i 



where K and i? are the complete elliptic integrals of the first and second kind respectively. The rest of the unit-order, 
odd-half-integer degree Legendre functions of the second kind can be computed using the following recurrence relation 

Notice that all odd-half-integer degree, integer order Legendre functions can be written as linear combinations of 
elliptic integrals of the first and second kind. The analogous formulas for q = are given in Cohl and Tohline (1999) 
[H], eqs. (22)-(26). 

Let us look at its behaviour for fi being a negative integer such that the binomial expansion should reduce to a 
polynomial in z. Let us take the limit as fi —q. Using the negative-order condition for Legendre functions of the 
second kind ([12, p. 367) 

, _ T{p-n) n+1/2, , 

and 

(, _ cos^)^ = „,y|(-i)«(,2 _ i),/2+i/4 g ,osH,)^^Qlt\%{z), (16) 
where we have used eq. ([6]). Let us explicitly verify 

Lq = lim(z — cos'0)'' = 1. 

q^Q 

Evaluating the limit of the right-hand side of eq. yields 

Lo = V ^»cosM) Tjn^ 

" V TT i ^ n\ ^"-1/2^ r(-g) ' 
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where we have used 



lim 3^ ^ 



Hence for p — 0, only the first term in the sum survives. The Legendre function 



is given in Abramowitz & Stegun (1964) (P, p. 334), reduces the Umit to 1, as expected. Take for instance 

Li = hm (z — cos tp)'^ — z — cos ip- 
Evaluating the limit of the right-hand side of eq. p6|) yields 

T /^"v 2 1^3/4X^ £nCOs(nVj) 3/2 ^{n-q) 

ji=0 ' ^ 



r(-9) 



The limit in the sum is given by 1 for rt = 0, — 1 for n = 1 and zero otherwise, hence only the first two terms in the 
sum survive. Using the Legendre functions 



1 /¥ 



-3/4 



we arrive at the desired result for Li. 

Substituting ^ = e N in eq. yields 



1 



cos l/']'^ 



'V n ^^Tv ^ cos(n^)g;'^_//2(z). 



(17) 



which the right-hand side is verified to be positive real upon examination of the definition of the Legendre function 
of the second kind in term of the Gauss hypergeometric function, eq. (|12p . For q = 1 the coefficients are Legendre 
functions with odd-half-integer degree and ^ order. These can be evaluated using ([T], eq. (8.6.10)) 



Ql"{z)^^^l{z^-l)-''' [z+Vz^l 
valid for S C and \z\ > 1. Note we also have ([1], eq. (8.6.11)) 



-1^-1/2 



If we take z = cosh rj and = n — i , where n E 1^ and insert the resulting expression in eq. (|17p we obtain 



1 



—— ^e„e-™'cos(n^). 



cosh rj — cos ■i/' sinh rj 
Similarly if we use eq. p9|) and the order recurrence relation for Legendre functions 



(18) 



(19) 



(20) 



)(;+2(z) = -2(/i + i)- 



Q^^+\z) + {ly-^^)(u + ^, + l)Q'^{z), 



we are able to compute all required odd- half- integer order Legendre functions appearing in cq. (I17p . 
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By taking z — cosh r\ we have 



1 



^ oo 

^ — 3 — e„ cos(nV')e^"''(nsinh?7 + coshr/), 



and 



(cosh 77 - cos -0)2 sinh^ 77 ^ 

1 ^00 

7 ; 7TT = cos(riV')e""'' 

(cosh 7/ - cos 2 sinh'' 77 ^ 

X ((n^ — 1) sinh^ 7; + Sti sinh 7y cosh 7; + 3 cosh^ 77) , 



00 

V e„ cos(7iV')e""'' 



(cosh 77 - cos 6 sinh" r\ 

X {[rC" — ^n) sinh"^ 77 + (671^ — 9) sinh^ 7; cosh 77 + 15nsinh77Cosh^ 77 + IScosh'^ 7;) . 

One way to verify these formulae is to start with the the generating function for Chebyshev polynomials of the 
first kind r„(a;) ([27] and [18], p.51) 

1 + Z2 _ 2X2 " ^ ^"^"^"(2^)' 
n— 

and substitute z = cosh 77, yielding eq. (j20p . The rest of the examples can be verified by direct repeated differentiation 
with respect to r/. 

3 Closed-form expressions for certain Legendre functions 

We have 

. / N p ( f^ + n fi + n+l _2 

An n[z) — ; \ — 2^1 , ; 71 + 1: Z 

' 2"7i!z'^+" V 2 2 ' ' ' 

-^^ ^rSw- (21) 

Let ^ 

2 = cosh (77) = - (x + a;""^) , 

with 

Since 77 > 0, we have < x < 1, and 

A. n (cosh 77) ^ MlL2'^X^+" (1 + x2)-(^+") 
TV. 

'/i + 71/i + 71+l 4x 

~' 2 + (TT"x2) 

Using Andrews, Askey, & Roy (1999) ( 2^ eq. (3.1.9)) with a = f.i + n and b = ^, we get 

A^,„ (cosh 77) - M« 2M^M+n 2^^^ (^ + n, /i; 71 + 1; x^) . 

71! ^ ' 

Using Pfaff's transformation ([2j, eq. (2.2.6)), we obtain 

A^,„ (cosh 7,) = Mil 2^x'^+" 2^^i (^i, + 71; n + 1; x^) 

^I^2^x^+"(l-x2)^%^^i ('/.,l-7i;n+l; ' 



X2^r(^,^^^^;n + l;,^). (22) 



C2 - 1 
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Taking /i = g G N, we can write 



or 



Ag^n (cosh 77) = 2« 



(1-X2)« 



9-1 

E 



fc=0 



q + n— l\fq + k — 1 



k 



9-1 

E 



V 2a; / fe=0 ^ 



a;2 - 1 
fe 



Since x = e ^, we get 



Aq^n (cosh 77) = 



-nr} 



sinh^ {rf) 



9-1 

E 

fc=0 



q + n— l\(q + k-l\ e 
71 + fc 



2*^ sinh 77 



or for instance by using eq.((TO] 



^'■"^'°'^''^-sinh^(7y)^^U-fc-Vl 9-1 y2^sinh^ 
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Using this formula and eq. (j2ip . we are able to write down a concise formula for odd- half- integer degree, odd-half-integer 
order Legendre functions 



2 (q- l)!(z2- 1)1/4 



9-1 

xE 

fc=0 

or in terms of Pochhammer symbols 



n + q — 1\ fq + k ~ \ 



q — k — 1 J \ (?— 1 



2(^2 _ 1)1/2 



yn-l/2l^)-H V2^^^ (z2 „ 1)1/4 



9-1 

E 

fc=0 



(g)fc(l - q)k 

{n+l)kk\ 



~z + \Jz^ - 1 

2(^2. 1)1/2 



Since 



Wln{q)k{^ -q)k _ fq + k-l\fq + n-l 

n + k 



{n + l)kkhi\ 

This leads us to conjecture the following general formula, q E N 

^ T{q + u+l) (z + ^/i^)--V2 



Qri/2(^) = ,(_1)9 



2 r(7/+|) (Z2_ 1)1/4 



9-1 

E 



{q)kil - q)k 



2(z2 _ 1)1/2 



This formula is verified using entry. (31) on p. 162 of 27J, g G Z to obtain 

/Fr((? + 7/ + i) (z-f Vi^^)— 1/2 



gri/2(^) ^ ,(_i)9^ 



2 r(i/-f|) (z2_ 1)1/4 



E 

/c=0 



(9)fc(l - 9)fe 



3^ 

2;fc ' 



-z -h Vz^ - 1 

2(z2- 1)1/2 



(23) 
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This is a generalization of the very important general formulae given by eqs. (8.6.10) and (8.6.11) in [T], p. 334. Taking 
z = cosh r] we have 

{q)k{l-q)k{-l)''e~''^ 



X 



Kq)k{-L - q)k {--L) c ■ 

^„ (i^+^)^kl 2fesinhS' 



k=0 



An alternative procedure for computing these Legendre functions is to start with eqs. (|19p and (|18p and use the order 
recurrence relation, [52] eq. (8.732.3) 

— 1 

On the other hand, the expressions given by eqs. (j23p and (|24p directly give closed-form expressions for the Legendre 
functions simply by evaluating a finite sum. 



4 Conclusion 

Not only is the generalized Heine identity useful for studying Poisson's equation in three-dimensions, but it is equally 
valid with fundamental solutions for k powers of the Laplacian in M". These are given in terms of a functional form 
which matches exactly the generalized Heine identity, particularly those in odd dimensions and in the even dimensions 
for k < ^ — 1. As will be seen in follow-up publications, the generalized Heine identity can be used as a powerful 
tool for expressing geometric properties (multi-summation addition theorems) for these fundamental solutions in 
rotationally-invariant coordinate systems which yield solutions to such equations through separation of variables. 
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